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In a previous paper of one of us [Europhys. Lett. 59 (2002), 330-336] the validity of Greene's 
method for determining the critical constant of the standard map (SM) was questioned on the basis 
of some numerical findings. Here we come back to that analysis and we provide an interpretation of 
the numerical results by showing that no contradiction is found with respect to Greene's method. 
We show that the previous results based on the expansion in Lindstedt series do correspond to the 
transition value but for a different map: the semi-standard map (SSM). Moreover, we study the 
expansion obtained from the SM and SSM by suppressing the small divisors. The first case turns 
out to be related to Kepler's equation after a proper transformation of variables. In both cases 
we give an analytical solution for the radius of convergence, that represents the singularity in the 
complex plane closest to the origin. Also here, the radius of convergence of the SM's analogue turns 
out to be lower than the one of the SSM. However, despite the absence of small denominators these 
two radii are lower than the ones of the true maps for golden mean winding numbers. Finally, the 
analyticity domain and, in particular, the critical constant for the two maps without small divisors 
are studied analytically and numerically. The analyticity domain appears to be an perfect circle for 
the SSM analogue, while it is stretched along the real axis for the SM analogue yielding a critical 
constant that is larger than its radius of convergence. 

PACS numbers: 05.45.-a, 05. 45. Ac, 45.10.Hj 



I. INTRODUCTION 

The Taylor-Chirikov map Q, @ or standard map (SM) 
is one of the best known nonlinear models showing the 
onset of chaos in Hamiltonian systems. It describes 
with some level of approximation many physical systems. 
Among these there are numerous applications to plasma- 
physics, in which field it was originally introduced. The 
SM is also exactly related to the time evolution of the 
"kicked rotor" and the equilibrium condition for a chain 
of masses superpositioned on a periodic potential. The 
latter is known as the Frenkel-Kontorova (FK) model. 
This model is of equally importance for solid state physics 
as the SM is for plasma physics. It has, e.g., been 
applied to Joscphson junctions arrays, charge density 
waves and surface friction 3]. More importantly, due 
to the simplicity and, yet, their complex behavior, these 
minimalistic models have had an enormous impact for 
our understanding in complex phenomena such as non- 
linearity, chaos, quasi-periodicity, and commensurate- 
incommensurate transitions. Although now part of any 
text-book in nonlinear physics and studied extensively 
over many years, the SM and FK still bear many un- 
solved problems. The most intriguing one of these is the 
sudden transition from smooth to chaotic orbits in the 
SM when the coupling parameter K is increased above 
a critical value K c . In the FK model this transition is 
connected to change from a sliding to a pinned state and 
bears the name Aubry- or breaking analyticity- transi- 
tion. 

The theoretical framework that characterizes this 



transition originates from the Kolmogorov-Arnold-Moser 
(KAM) theorem Q that deals with the problem of small 
denominators that can occur in any perturbation expan- 
sion for quasi-integrable systems. In fact, the KAM the- 
orem can be used to prove the non-chaotic behavior of 
the SM for very small coupling K and sufficient irra- 
tional winding number I. Other arguments can then be 
applied to prove that a chaotic regime exists for values 
K > K' giving an upper bound to K c . For I equal to 
the golden mean, there exists an analytical bound by 
Mather K c < 4/3 [jj, and the computer assisted proof of 
MacKay and Percival K c < 63/64 w 0.9844 @. More- 
over, another computer assisted analysis of Jungreis ex- 
cluded the value K = 0.9718 forpossible occurrence of 
invariant circles (smooth orbits) [3. 

There exist several methods to calculate K c precisely, 
among which Greene's method Q has shown to be one 
of the most effective giving the estimate K c = 0.971635. 
This method is based on the assumption that the disso- 
lution of invariant curves can be associated with the sud- 
den change from stability to instability of nearly closed 
orbits. The renormalization method of MacKay Q is 
a further refinement of this method and has established 
the same value with higher digit precision with respect 
to the original Greene's result. Yet, Greene's hypoth- 
esis has only been partly proven. A result by Falcol- 
ini and de la Llave |9( asserts that the critical constants 
for symplectic maps can never be higher than the ones 
obtained by Greene's method. Recently, the result has 
been extended to nontwist maps by Delshams and de la 
Llave :ifj|. Hence, K c < 0.971635 for the SM with the 
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golden mean as winding number. 

Another way to calculate K c is through the Lindst- 
edt series expansion. The smooth orbits in the SM can 
be described by a continuous (analytic) periodic function 
beneath K c . Hence, below K c the Fourier spectrum of 
this function should be finite for each component while 
decaying to zero in the high frequency domain. Above 
K c some divergence is expected: either infinitely high 
frequency components persist, or the amplitude of some 
components in the spectrum diverges, or both. By writ- 
ing down the Taylor expansion and equating the Taylor 
orders in the functional equation satisfied by this conju- 
gation function, the Fourier- Taylor components can, in 
principle, be derived from the ones of lower order. In 
Ref. an evaluation of this expansion to high orders 
suggested a convergence to a value K c ~ 0.97978, which 
is higher than Greene's result. In this article, we come 
back to this analysis and show that an apparent plausible 
assumption made in Ref. is falsified beyond Taylor 
order n > 200. As a result, the Lindstedt expansion does 
not contradict Greene's result. The value K c ~ 0.97978, 
however, does correspond to the critical value for a differ- 
ent map, the semi-standard map (SSM). We come back 
to this in Sec. IIIII 

Aubry ^| proposed another method, that is probably 
not very effective for high precision evaluation in a com- 
puter algorithm, but still interesting. It is based on an 
eigenvalue calculation of the dynamical matrix for the FK 
chain close to the critical point. Although this, in princi- 
ple, requires the diagonalization of an infinite matrix, one 
can use the fact that the eigenvector of the lowest mode 
tends to localize 0] . The instability of the FK chain can 
than be determined in successive approximants by calcu- 
lating the determinants of finite matrices of increasing 
length. 

Finally we mention the use of Pade approximants |3Sj 
to study numerically the entire analyticity domain. This 
is a powerful tool even if is it less precise than other 
methods for detecting the critical constant K c and not 
completely under control from a rigorous poin t of view. 
It has, for instance, been employed in Ref. [lj| and, very 
recently, in Ref. where the existence of a natural 
boundary for the analyticity domain of the SM has been 
checked numerically Always with the aim of st udy ing the 
analyticity domain Falcolini and de la Llave [lg] devel- 
oped a variant of Greene's method working for complex 
values of the parameter K that gives an alternative to 
the Pade series approach. 

Eventually, these approaches are assumed to converge 
to the same value. However, the proof of this is highly 
non-trivial. The ultimate goal, of course, would be to 
gain an analytical expression for K c . This is still far 
beyond our capabilities. Inspired to investigate further 
the influence of the small denominators in the Lindstedt 
series expansion, we introduce two simplified models by 
setting rigorously all the denominators equal to one both 
for the SSM and the SM. In the latter case, this is a very 
well known model, Kepler's equation |l7j . which turns 



out to have a very similar transition and can be solved 
analytically. The radius of convergence is found to be 
higher than the SM and SSM value in case of golden 
mean winding numbers. 

This article is organized as follows: In Sec. [H] we intro- 
duce the SM and SSM. In Sec. 11111 we come back to the 
analysis of Ref. ^l] showing that the Lindstedt expansion 
does not violate Greene's method and make the compar- 
ison between the SM and SSM. In Sec. IIVI we present 
a new model in which we suppress the small denomina- 
tors and give an analytical expression both for the radius 
of convergence and the critical constant. We end with 
conclusions in Sec. 



II. THE (SEMI-) STANDARD MAP 

The SM and SSM are defined as 

Xi \ ( 2xi + V'(xi) - Xi-i 



Xi+l 



T 



Xi-l 



with 



V'{x) =^sin(27nr) 
v '( x ) = as exp^Tra) 



for the SM, 
for the SSM. 



(1) 



(2) 



The resulting sequence (xi mod 1) for i = 2, . . . oo origi- 
nating from a starting point (x\, xq) produces a discrete 
trajectory on a two-dimensional torus. Such a trajectory 
for the SM can be related to the equilibrium positions 
of an infinite FK chain where particles with harmonic 
nearest neighbor coupling are placed on a periodic po- 
tential V{x) = K{2ir)- 2 {1 - cos(2ttx)). The SSM has 
not such a similar counterpart, but is much simpler in its 
mathematics. By definition, I = (Xi+i — Xi) is called the 
winding number or rotation number of the map. For low 
coupling K and I incommensurate to the periodicity of 
V , there exists a continuous function gi(x) such that the 
positions {xi} can be expressed as 



Xi = gi(il + a). 



(3) 



This function is often called the conjugating function or, 
in context with the FK model, the modulation or hull 
function. The subscript I indicates that the shape of the 
function depends on the rotation number. The conjugat- 
ing function satisfies the functional equation 

2gi{x) -gi{x + l)- gi(x - I) = -V'(x + gi(x)) . (4) 

For K large enough the function gi(x) becomes discon- 
tinuous. For the SM this implies that the orbits become 
chaotic and for the FK that the chain of particles gets 
pinned together with the appearance of a phonon gap. 
This transition, in context with the FK model, is also 
called analyticity breaking transition or Aubry transi- 
tion. 

There are several quantities of interest which one can 
introduce in order to study the transition from regular to 
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chaotic dynamics. As the function gi(x) is analytic for K 
close to the origin one can consider its series expansion 
in powers of K 



9l{x) 



and define the radius of convergence p{l) as 



(5) 



Besides only partly proven, Greene's method has also 
some other limitations. For instance, this method does 
not work for other interesting models, as the SSM and 
Siegel's problem |l8[ . where the construction of periodic 
orbits fails |4l| . The best general alternative is the Lind- 
stedt series expansion. This method is more generally 
applicable (it also works for the SSM), but is less accu- 
rate than Greene's method for the SM. 



p(l)= inf (hmsup| 3 r(z)| 1/Tl ) \ (6) 

xe[0,l] V n-»oo / 

where infzefo.i] denotes the infimum or greatest lower 
bound in the domain x S [0, 1] for the quantity within 
the brackets (. . and limsup 1|woo is the supremum 

. , \l/k 

limit giving the highest value for \gf (x)\ of all k > n 
in the limit n — ► 00. Note that the infimum appears in 
the definition of the radius of convergence because each 
invariant curve is filled densely by any trajectory lying 
on it as a result of the incommensurate winding number. 
Hence, existence of the invariant curve itself requires the 
latter to be defined for all x E [0, 1]. 

The critical constant K c (l) is defined as the (positive) 
real value K c (l) such that for K > KM) the conjugat- 
ing function is not analytic any more 39] . It is believed 
that the analyticity domain of the conjugating function 
has a natural boundary: this means that g(K,x), which 
is the modulation function g{x) at a value K , has a set 
of singularities in terms of K that form a closed curve 
around the origin in the complex plane. Hence, the ra- 
dius of convergence p corresponds to the singularity clos- 
est to the origin, while the critical constant K c corre- 
sponds to the intersection of this curve with the (posi- 
tive) real axis. By definition one has K c (l) > p(l), so 
that by estimating the radius of convergence one finds 
a lower bound for the critical constant. Furthermore it 
is generally accepted that K c (t) = p(r) for the golden 
mean r = (\/5 - l)/2 w 0.618034 0. It is also com- 
monly believed (on the basis of numerical simulations and 
heuristic arguments) that K C {1) has the highest value for 
the golden mean r. Moreover, KM) is assumed the same 
for all values I S Z(r), i.e. the values that can be written 
as I = mr + n with m, n integer numbers. 

So far, the most accurate method to calculate K c is 
based on Greene's method (also known as residue crite- 
rion). In this method the infinite trajectory {x{\ with 
irrational winding number I is approached by successive 
approximants j, which are periodic trajectories with ra- 
tional winding numbers lj — pj / qj and Xi+ q . mod 1 = x±. 
Hence, qj and pj are at each level j two integer values 
giving a better estimate of I for each increment in j and 
I = limj—xjo lj. These numbers can for instance be ob- 
tained using the continued fraction fraction expansion 
of I. For I = t this results in the Fibonacci numbers 
(r « Fj-i/Fj with F = F x = 1 and Fj = Fj-i + Fj for 
j > 1). Conclusively, the Greene's method tells how to 
construct the periodic orbits and to measure their stabil- 
ity by means of the residue that does not tend to zero 
any more for K > K c . 



III. THE LINDSTEDT SERIES EXPANSION 

Standard Map: A way to study the transition is by 
means of the Lindstedt series, which is the expansion of 
the function gi(x) both in Fourier and in Taylor series 
[42j. By defining the Fourier transform as 



9l{x) 



X k = 



X k e 27rtkx with inverse: 



dx gi(x)e 



and expanding 



X k = KXl + K 2 Xl + K s Xl 



(7) 



(8) 



we end up with Fourier- Taylor coefficients X2, where n is 
the Taylor index and k is the Fourier index. Now, using 
Eq. we can relate the Fourier- Taylor coefficients of 
order n by the ones with lower Taylor index by |ll| 



1 1 2 Y n 



-\5l,k - <S-l,fc}<5l,n H 



47T I 



47T 



^ (i27T) m 

/ T~ x 



E { E 

ni+U2 + ...+n m — n— 1 /ci+A;2 + ---+A; 7T1 — k — 1 

-(-ir E 
/ci+fc 2 +...+fc m =fc+i 



A fci A k 2 



A fel X k 2 



■ XI 



with 



Xk Jo 



dx (2g l (x) -g t (x + l)- g t (x - I)) e 



= 2(1 - cos(27rfc0) = (2sin(7rfc0)' 



-2-nikx 



(10) 



and where V" , . , , implies a summation of 

all possible integers ni, n^, ■ ■ ■ , n m with the constraint 
that 53j=i n i = n — 1. There are ways to reduce the 
number of summations in Eq. (|5J). One possible way was 
proposed in Ref. jl to construct an extended matrix 
P(n, k, m) defined as 



P(n, fc, to) = 

E 



{2m) 7 



E 



ni +n2 +■ ■ ■ -\-n m —n 



ki /C2 



■ XI 



(ii) 



k 1 +k 2 + ...+k m =k 
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One can show that P(n, k,m) =0 if \k\ > n or m > n. 
This gives rise to following recursive relations (lTI| : 



Fi-i 



F 3 r 



{— t) j ' +1 , one can show that for j — > oo 



P(l,±l,l) 
P(n,fc,l) 



2w 2 ' 



(12) 



n — 1 

-^X^ [-P(n-l,fe-l,m) 

rn — 1 

(-l) m P(n- 1,*+ 1,171)1, 



n— rn+1 min{n',fc+n— n'} 

P(n,fc,m) = - £ ^ 

n' — 1 fc'-maxj-n' n+n'} 

P(n', fc', l)P(n -n', k- k', m - 1), m > 1, 

from which we can distract the Fourier- Taylor coefficients 
by 



P(n,k,l) 
2m 



(13) 



The components of the P matrix are all real and obey 
the symmetry relation P(n,k,m) = (—l) m P(n,—k,m). 
Moreover, besides being zero for \k\ > n and m > n, 
the matrix P(n,k, m) has zero values whenever k + n is 
odd. Hence, k = n, n — 2, . . . , —n are the only nonzero 
elements of P. 

The evaluation of Eqs. I|12fl is very efficient to evaluate 
X% obtaining a Taylor order of approximately n = 200. 
To go beyond this limit, sufficient computer power and 
time is needed as both the computation time as the 
number of nonzero matrix elements increase with ~ n 3 . 
Hence, memory can become a severe problem as the num- 
ber of components that have to be stored can easily go 
beyond the maximum allowed allocation limit of the com- 
puting system. In this work, we finally reached the level 
n = 700 (see Fig.^l and we believe that going beyond this 
order is not very profitable for obtaining a more accurate 
evaluation of K c . We come back to these results after ad- 
dressing the small denominator problem that arises from 
Eq. 0. 

From Eq. (|10|l and the second line in Eq. I|12|) one 

sees that even for irrational values of I, the terms uj^ 2 
can become arbitrarily high for some k. This effect is a 
typical example of the 'small denominator problem' that 
can strongly prevent the convergence of any perturbative 
series. In fact, it requires a stronger condition than ir- 
rationality such as a Diophantine condition Among 
all the irrational numbers, the golden mean r suffers the 
least from the small denominator problem and has there- 
fore the highest convergence radius K c . The golden mean 
is relatively difficult to approximate by rational numbers 
as those arising, for instance, from the continued fraction 
expansion (one can say that it is the 'most irrational 
number'). As addressed above, the best approximants 
for the golden mean are given by the Fibonacci numbers: 
r = Fj / Fj + \. Therefore, the small denominators will 
be most dominant for k = P. . From the exact relation 



u> 



4tt 2 (t 2 ) 



3+1 



(2.618) 



j'+i 



(14) 



However, we would like to stress that the small denom- 
inator problem is not the only mechanism causing the 
breakdown of the perturbative approach. This becomes 
evident in Sec. IIVI where we introduce the model that 
arises when we rigorously set lu^ 2 = 1 for all k in the 
series of Eq. I|12|) . Clearly, this simplified expansion can 
not be affected by the small denominators. However, it 
still has a radius of convergence and a critical constant, 
as shown by the analytical solution. As the radius of con- 
vergence p of this simplified model is found to be lower 
than p(r) — K c (t) for the SM, it proves that the golden 
mean winding numbers I G Z(t) are remarkably resistant 
to the problem of small denominators. The full analysis 
of this model is given in Sec. IIVI 

Coming back to the results of Fig.^ we see that indeed 
the evolution of P(n, k, 1) makes sudden jumps at the 
Fibonacci numbers as expected from Eq. (| 1 41) . Besides, 
deceptively the values n = k = F n seems to determine 
the whole power law behavior of XJ}, which is true till 
Taylor order n ~ 200. This assumption made in Ref. pl| 
allows for a further simplification of Eq. I|12|) by defining 
the reduced matrix Q(n,m) = P(n,n,m) obeying the 
relations 



Q(i,i) 

Q(n,l) 
Q{n, m) 



-1 
2^' 



-^- 2 J2Q(n-hm), (15) 

in— 1 
n— rn+l 

— Q(n',l)Q(n-n',m-l). 



This set of equations make high order (~ n = P20 = 
10946) evaluations accessible for computer calculations. 
At this order the value of the radius of convergence seems 
to stabilize near p = 0.97978 , but this is still higher than 
the one obtained by Greene's method. As consequence, 
validity of Greene's method was questioned in the quoted 
paper pT| . 

The more elaborated calculations in this work show 
that the assumption made in Ref. |ll| was actually wrong 
as shown by the high order behavior in Fig. ^ Still, we 
find that k — n gives the maximum for |P(n, fc, 1)| when- 
ever n is a Fibonacci number. However, the character of 
the evolution changes from being peaked to more smooth 
oscillations. Clearly, the line connecting \P(Fj, Fj, 1)| 
does no longer dominate the increment of P-matrix el- 
ements for n > 200. As |P(n, fc max (n), l)! 1 /™ shows lo- 
cal maxima for (n, fc max ) = (383, 377) and (n, fc max ) = 
(622, 610) just after P13 and P14, we fitted the line o^A^ 
through these points. From this fit, A2 = 1.0248, the 
estimate for K c « A^ 1 = 0.9758 is obtained. Although 
still higher than Greene's value, it is already consider- 
ably lower than p — 0.97978 obtained from Eqs. i|15[l for 
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n = F 2 n = 10946 0] • Note that the latter approach of 
Ref. [llj for this lower approximant n = F 14 = 610, as 
obtained from the line A" (See Fig. [JJ, would result in 
psi 1/Ai = 0.9817, still approximately 0.002 higher than 
the nearly converged value of p = 0.97978. Hence, a de- 
cay of 0.004 from 0.9758 at n = 622 to 0.9716 at n -> oo 
is not unlikely. As a consequence, contrary to the results 
of the restricted series 1)150 , there is no proof that the full 
Lindstedt series (|12ll violates Greene's hypothesis. This 
also shows that further simplification of Eqs. II 1211 is not 
easily obtained and that an accurate evaluation of K c 
based on the Lindstedt perturbation is severely demand- 
ing. 




n 



FIG. 1: \P(n, fc max , 1)1 as a function of n. This is denned 
as the maximum value of \P(n,k, 1)| of all k. Hence, fc max 
is defined as the k value where \P(n,k,l)\ has this maxi- 
mum. The inset in the lower corner shows fc max as function 
of n. The inset in the left upper corner is an enlargement 
of the first 200 terms together with \P(n,n, 1)| (dashed line). 
From these figures, one can clearly detect sudden boosts in 
the function \P(n, fc max , 1)| where fc max = n at the Fibonacci 
values (dashed vertical lines). However, whereas for n < 200 
the character is sharply peaked at these values, its behavior 
changes for higher orders. Still fc max = n for n a Fibonacci 
number, but the intersecting line described by aiA™ does no 
longer dominate the complete evolution of all the \P{n, k, 1)| 
terms. Ai = 1.0186 is determined by the line through (n, k) — 
(Fi 3 ,Fi 3 ) = (377,377) and (n, k) = (F 14 ,F 14 ) = (610,610). 
A 2 = 1.0248 is set by the line through (n,k) = (383,377) 
and (n,k) = (622,610) where \P(n, fc max (n), l)| 1/n shows lo- 
cal maxima in n. The inversed values, A^ 1 « 0.9817 and 
Aj 1 « 0.9758 are assumed to converge for higher n to K c for 
the SSM and SM, respectively. 



Semi-Standard Map: When evaluating Eq. Q for 
the SSM ©, the factors <J_ 1)fc and -(-l) m E • • • are not 
present in Eq. It is then straightforward to show that 
the P- matrix Ijlll) has only non-zero elements for n = k. 
Hence, the assumption made in Ref. that gave rise 
to Eqs. I|15|) does not corresponds to the critical constant 
of the SM, but still gives the correct value for SSM. A 
result by Davie [n| shows that, for maps including the 
ones we are considering, the radius of convergence © is 



equal to 

p(l) = (limsupmaxlA^I 17 ™) . (16) 

\ n— >oo \k\<n ' 

Therefore, the radius of convergence for the SM cannot 
be larger than the radius of convergence of the SSM, but 
of course it implies only a lower bound on the critical 
constant. Numerically by using Pade approximants in 
Ref. 0] it has been found that for certain values of the 
rotation number I, the radius of convergence of the SM 
is strictly smaller than the radius for the SSM. For the 
golden mean it is hard to improve upon simple power se- 
ries using Pade, but for other numbers closer to resonant 
values it is possible and the phenomenon becomes much 
more evident. 

The fact that the critical constant for these numbers is 
lower for the SM than for the SSM implies, by fity. that 
dominant contributions arise from terms with Taylor or- 
ders n for which |X£ | > |A™|. This is exactly what 
emerges from the numerics as noted above and shown in 
Fig-fflfor n > 200. Clearly, this is not be the case for the 
SSM where one can limit to k = n. Hence, in Ref. 0] 
K c {t) = p{r) was actually determined for the SSM to be 
0.97977 at Taylor order n = F 20 = 10946. In the calcu- 
lations of this work, we went to order n — F24 — 75025, 
that gave the value 0.97937. As the root criterion sat- 
urates very slowly, the numerical results provide essen- 
tially only an (accurate) upper bound for the radius of 
convergence. 

To conclude this section, we found that, also for the 
golden mean, the radius of convergence for the SM is 
strictly less than the radius of convergence for the SSM. 
Therefore, as a general comment we can remark that for 
the SM the presence of all harmonics in the Fourier ex- 
pansion of the coefficients gV(x) has a double effect. On 
one hand the radius of convergence becomes smaller with 
respect to that of the SSM. On the other hand, the crit- 
ical constant K C {1) can be larger than its radius of con- 
vergence p{f) for I £ Z(t). For the golden mean the 
two values are equal as emerges numerically , but for 
other values they can be appreciably different. One can 
imagine that the first phenomenon is due to the pres- 
ence of contributions larger than X™, while the sec- 
ond one is a consequence of deep cancellations between 
the harmonics of given perturbative order. These two 
effects are, in general, much more dominant for rota- 
tion numbers I close to rational values (see for instance 
Refs. [Hill El). 

IV. PUTTING SMALL DENOMINATORS TO 
UNITY 

A. Introduction of the simplified maps 

An interesting study appears if we set rigorously all 
possible small denominators equal to unity: u>k = 1 for 
all k in Eqs. I|12fl and Q15[l. Although the inspiration of 
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this model was simply to study of the perturbation ex- 
pansion when the small denominators have no effect, we 
can retrace from this series back to a functional relation 
as the one in Eq. ft can be shown that these sim- 
plified series would correspond to the functional relation 
for a function h(x) 

K 

h(x) H exp (i2Tv(x + h(x))) = for the SSM and 

h(x) + — sin (2vr(a; + h{x)) = for the SM. (17) 

Hence, the divergence of the simplified series corresponds 
to the point in K where these functional forms (|17|) have 
no analytical solution any more. A logical next step 
would be to relate the equality (|17fl to the iteration of 
a map similar to Eq. Q. As the relation (|17l) does no 
longer contain the arguments x± I this is not so evident. 
However, one can relate the h(x) function to the hull 
function of a FK-type system. It can be shown that this 
corresponds to an one-dimensional Einstein solid that is 
interacting with an external incommensurate potential. 
Due to the lack of neighbor interaction, which makes each 
particle independent, it is highly unusual to describe for 
such a system the equilibrium coordinates by a collective 
hull function. Still, there are no restriction not to do so 
and one can even give such a function a physical meaning. 
As known from the FK model, the continuous shape of 
the hull function is directly associated with the existence 
of a sliding mode where the FK chain can slide over the 
periodic potential without cost of energy 0, ^| . In this 
case, the complete phonon spectrum is giving by sum of 
oscillations of the individual particles that are not zero 
in general. The sliding mode appears when we add an 
extra degree of freedom to the system as shown in Fig. EI 




< 5 




FIG. 2: Illustration of the FK model (top) and the system 
that obeys Eq. 1171 (bottom). The latter corresponds also to 
the FKT model without neighbor interaction. All particle are 
connected to the upper rod whose position is given by £. A 
sliding mode may exist when £ can be varied. 

Here, all particles have no interaction to their neigh- 
bors, but are connected to an upper rod. When the rod 
has an infinite mass compared to the particle masses, the 
system is basically an Einstein solid. However, if we as- 
sume that the position of the rod may vary according to 



a coordinate £ an extra phonon mode exist that is zero 
for K < K c in this system. Hence, the breakdown of the 
Lindstedt expansion (|12fl with Uk — 1 for all k can be also 
related to a real physical sliding-pinning transition. The 
model as illustrated in Fig. is also equal to a special 
case of the Frenkel-Kontorova-Tomlinson (FKT) model 
where each particle is connected to a rod by a spring with 
spring constant c v and additionally to its neighbors with 
a coupling c„ |43- The FKT model has been proposed 
to study more realisticl y th e frictional behavior between 
atomic surfaces |22ll23.l2^| . The model that is described 
by Eq. (|17l) simply corresponds to the FKT system with 
c r = 1 and c„ = 0. The nice thing is that the perturba- 
tion series of Eqs. (|T7|) can be solved exactly. To show 
this, we will start with the more simple SSM case. 

B. Radii of convergence for simplified maps 

Simplified SSM: It is convenient to use following nor- 
malization 

R(n,m) = n\(-2) n Q(n,m), (18) 

with matrix entries that are integer and positive and with 
R(n,n) = 1. From Eq. (JT2J| and JTSJ with lu~ 2 = 1 we 
derive 

ra-l 

R{n,l) = n^2R(n-l,m), (19) 

7/7 = 1 

R(n,m) = — V )R(n',l)R(n-n',m- 1), 

n' = l v ' 

with ( n .) the binomial coefficient: = ,,, " ! — ttt. Note that 
the recursive relations in Eq. i|19|) f° r m > 1 coincide 
with those satisfied by the Stirling numbers of first and 
second kind, Sn and 6^ respectively Of course 
what is different is the relation for m = 1. 

From relations (|19|l with R(l, 1) = 1 the following ex- 
act equality can be proven 

R(n,m) = n n -" l ( U ~ 1 \ (20) 
\m — 1/ 

The proof of this equation is given in the appendix in 
two ways. In App. IA II we derive this proof using the 
argument of induction. In App. IB II we give a proof 
based on the tree formalism that was firstly introduced 
in Ref . [25| . The first proof is straightforward, but elon- 
gated. The second is short, but less self-contained as 
it requires some knowledge of the previous publications 
about the tree formalism. The latter is, however, most 
practical for the more complicated proof of relation for 
the SM in series of Eq. I)17|l with lo^ 2 = 1. Now, from 
Eq. and J50I) we deduce that 

S(n,l) = (-ir|^, (21) 
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which is assumed to increases as a power law ~ a\ n giv- 
ing the radius of convergence as p = 1/A. Hence, 

In | Q(n, 1)| = (n- l)ln(n) - nln(2) - ln(n!). (22) 

Then using (2(| 

ln(ra!) = (n+^jlnn- n+^\n(2n), (23) 

that is an refinement of the well known Stirling's formula 
In n\ « n In n — n, we get 



ln|Q(n,l)| = n(l-ln(2))-^ln(n)-iln(27r) 



|Q(n,l)| 



1 (\ 



V27m 3 \2 



(24) 



yielding a radius of convergence p — 2/e w 0.735759. 
This value is less than the SSM value or SM value. This 
is a bit contra-intuitive, as one would expect that the 
possible occurrence of small denominators would give a 
lower p. Apparently, this does not happen for the golden 
mean. This can be understood by following reasoning. 
Although the small denominator factors cu^ 2 can become 
arbitrary large for some k giving a boost to the series 1(1 2|) 
and (|15J) . at most values of k they will be considerable 
smaller than 1 resulting in an opposite effect. Hence, 
for the golden mean as winding number the latter effect 
seems to be more dominant yielding an even higher value 
for p than the case where all w^T 2 terms are equal to 1. 

Simplified SM: The simplified SM considered in 
Eq. 1(17|) is well known in celestial mechanics 0] after 
applying the following variable transformation. Write 
K/2 = — c and 2irx — M where e is the eccentricity and 
M is the mean anomaly. Then the eccentric anomaly 
E — 2tt(x + h(x)) is related to M through Kepler's equa- 
tion M = E— esinE, which is exactly the second equation 
in Eq. fTfl . 

The recursive relations 1(12(1 with lo^ 2 = 1 have also an 
exact solution that we write here: 



Eq. 125(1 with Stirling's formula 1(23(1 . This gives rise to 

1 un—l p n 

I P f /t 1 1 I 

~ 2™ (\(n - k))h^- k ){\(n + jfe))*( n + fc ) 

' (26) 



(I(n2-fc2))l 
1 



■ 4 

2e 



(l-cr^-^l + O-)!^) 



where we have defined a — k/n € [—1,1]. Hence, we 
have to compute the maximum of the function 



E(a) 



a exp 



1 — a , , \ 1 + 0",. 
.— -ln(l-a) — ln(l + o) 



By taking the derivative ^ 



(27) 

0, we find that the max- 
imum is reached at a value o max that satisfies following 
relation: 

2 + (Tmax ln(l - o max ) - er max ln(l + cr max ) = 0, (28) 

yielding cr max w 0.833557. Hence, /c max = (J m ^n w 
0.833557n. Using Eq. QBJ, S(cr max ) simplifies to 
£J(o max ) = 1— j=ss=. Inserting this relation into 



Eq. gives 



|P(n,fc„ 



vn+l 



(29) 



with A = 



This yields a radius of conver- 



gence p = A -1 ss 0.662743, which is known as the 
Laplace limit [2£|. This value is again smaller than the 
radius of convergence of the true SM (recall that for 
the golden mean the radius of convergence p(r) equals 
K c (t)). Moreover, similar to the true maps, this SM- 
analogue transition value p is lower than the one of the 
SSM. 



( _-\\n+(n—k)/2 un— 1 

[ ' ' ' 2" ((n-fc)/2)!((n + fc)/2)!' 

for|fc| < n and if k + n is even 
= otherwise, (25) 

which can be obtained by the Lagrange inversion theo- 
rem |0, [23) . We present a new derivation of this relation 
based on the tree formalism in the App. IB 21 

Then, by using Eq. i|16fl we see that we have to com- 
pute the maximum over k of \P(n, k, 1)|. By assuming 
that the maximum is reached for some k which is not too 
close to n (an assumption that we shall verify a posteri- 
ori |45|). we can approximate the factorial appearing in 



C. The critical constant and the analyticity domain 

The argument above gives only information about the 
location of the singularities closest to the origin. The 
solution of the functional equations (|17fl could still exists 
for real values of K larger than the radius p. This is the 
situation K > p where h{K,x) is still analytical, but at 
which the power series (JSJ and JSJ) are no more defined. 
This is typical for a summation that consists of both 
positive as negative terms. In particular, there could be 
no singularity at all on the real axis so that an analytical 
form of h(K, x) could still exists for K — > 00. 

To analyze the extent of the analyticity domain and the 
critical constants, we need to 'evaluate' the summations 
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of Eqs. 10 and Q. This means that we need to find 
functional form h(K,x) that corresponds to the power 
series, but, contrary to the summation itself, can still be 
perfectly defined for \K\ > p. 

In the following analysis, we will show that the ana- 
lyticity domains for the simplified maps are, like the for 
true maps, also constrained by a closed boundary. More 
precisely, we find that for fixed x there are only a few 
singularities, but the union over all x S [0, 1] of such 
singularities reconstruct a closed curve surrounding the 
origin. Hence, outside this natural boundary there is no 
function h(K, x) that can be obtained by an analytic con- 
tinuation of the power series around K = 0. Although 
very unlikely, this does not completely exclude the ex- 
istence of a very different function, say h(K,x), that is 
defined outside this domain and obeys Eq. (|17J) and may 
even persist for K — > oo. Recurrence phenomena of this 
kind are known to occur for certain maps |29j . but, for 
instance, this is not the case of the SM. 

Simplified SSM: As for this model, X k {K) = X%K n 
we can directly writing down the summation of Eq. (JSJ 
for h(K, x) function: 



h{K,x) = ^X k t K k e 2 * lkx . 



(30) 



fc=i 



Inserting the expression (|21[1 gives 

i 00 uk-i 



k=l 



1 oo 

— \Q(k,l)\K k e™ k ( 2x+ V (31) 



2ni 

k=l 

To find the full analyticity domain of h(K,x), one ba- 
sically has to fix a certain value for x, say x = x' , and 
search for the singularities in K of the function h(K, x') 
by e.g. using the Pade approximants method. Then, one 
has to repeat, in principle, this procedure for all possible 
values of x and collect the set of all singularities to con- 
struct the full analyticity domain. Finally, the radius of 
convergence p is then the complex singularity closest to 
the origin, while K c is the smallest (positive) real singu- 
larity, if any. 

Vice versa, we could also fix the argument <j> of the 
complex value K, such that K = \K\e 1 ^ . The summa- 
tion (|31|l will then be maximized for x = — (3^+5), where 
each term in the sum turns into a positive value. For 
these values of K and x, using the inclusion argument and 
the root criterion on the delimiting series, one can show 
that the radius of convergence is given by 2/e. Hence, for 
each x there is one singularity at K = _2 e _l27r!C ; and the 
complete set over all x forms a natural boundary that is 
a circle around the origin. 

Note that this is very different from the true maps. 
Although not proven, numerical studies (for instance 
Ref. [l5j and references quoted therein) suggest that for 
the SM and SSM the function g(K, x) has for each value 



of x, independent to its value, an infinite set of singulari- 
ties forming the same (for each x) natural boundary. Nu- 
merical analysis 15j shows that the natural boundary of 
the SSM is a circle, just as in this simplified model. This 
property appears to be true irrespective to the choice of 
I as long as it fulfills Diophantine condition Q. For the 
SM with golden mean as winding number this curve re- 
sembles close to a circle, but not very smooth and slightly 
elongated (about 1%) along the imaginary axis [l6| . 

Simplified SM: Taking the power series (JSJ for Xk 
for the simplified SM using Eq. (|25[) we have 



00 ^ 

MK) = J2—P(n,k,i)K 



2ui 

I (_l)«+(n-*0/ 2 



2™ 



7T- y k 

2m ^ 

n=|fc|,|fc|+2,... 

((n-fc)/2)!((n + *)/2)f 



Changing variables to j — (n — |fc|)/2 gives: 



(32) 



X k (K) 



(-1) 



±_ \ ^ 



2nik 2 2 J'+I fc lj!(j + |fc|)! 



\k\K) 



2j+\k\ 



(33) 



with J v (z) the Bessel function of the first kind |30l 131 
defined (for integers v) as 



(-1)1 



3=0 



(34) 



As these Bessel functions J v (z) have no singularities in z, 
neither has Xk(K) in K. Therefore, the Fourier compo- 
nents do not give direct information about K c . On the 
other hand, one can conclude from Eqs. (|33|l and H34|) 
that |Afe(_ftT)| is maximized for pure imaginary K, so that 
the radius of convergence is lying on the imaginary axis 
on a distance p from the origin. Here, the individual 
terms in Eq. 134fl can not cancel as (—1)-' is then neu- 
tralized by x 2j ~ K 2j = (— l) J \K\ 2 i . Hence, we expect to 
obtain the singularity equal to the radius of convergence 
for K along the imaginary axis. 

We can now try to evaluate the Fourier series J7J for 
h(x): 



+ OO 



h(K,x) = 



(~l) fe 
2nik 



Jifei (l&l-K") exp(27rifca;) 



k— — 00 



= V ~-^-J k (kK)sm(2irkx). (35) 

^ — ' Trie 



k=l 



nk 



Further simplification is achieved by taking the derivative 
to x and searching the singularities in 

00 

h'(K,x) = ^2(-l) k J k (kK)cos(2Tikx) (36) 
fc=i 
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instead of h(K,x). This is allowed as the two problems 
are equivalent. 

From the series (|36|) we can guess for which values of x 
the singularities will be K c and p respectively. As Jk(kK) 
is positive for real values < K < 1 (see p. 534 in |3lj . 
we need to compensate the (— l) fe term by cos(27rfca;). 
This is achieved for x = h that reduces Eq. (|36[) to 

1 °° 
h'(K,-) = 2Y,MkK), (37) 

fc=i 

which has the exact solution (see formula (1) at p. 615 
in H2) 

°° K 
2Y,Jk(kK) = T -^. (38) 

k=l 

Hence, h'(K, ~) has a singularity at K — 1, yielding the 
critical constant K c — 1, a well known result in celestial 
mechanics [lTj. 

The complete analyticity domain can be found in 
Ref. 0], P- 219. In Fig. we represent what can be 
obtained by using Pade approximants for some values of 
x. What emerges is that the function h(x,K) has for 
each value of £ a pair of complex singularities closest to 
the origin symmetric with respect to the real axis. For x 
going from to 1/2 such singularities move continuously 
from —1 to 1 along two (symmetric) curves which pass 
through the points ±ip at x = 1/4 (see Fig. EJ). Hence 
the entire set of singularities closest to the origin lies on 
a curve which is smooth except at K = ±1, where it has 
a discontinuity in its first derivative (cf. again Ref. |l7j . 
p. 219). An important feature is, however, that, as al- 
ready noted in a similar context by Simon 32], a natural 
boundary in K for fixed x seems to appear only in the 
presence of small divisors. In fact, the latter give rise to 
the occurrence of sudden peaks yielding a pattern simi- 
lar to lacunary series |33j | for which natural boundaries 
can be proved to arise. Hence, these peaks seem to be 
responsible of the formation of the natural boundary as 
suggested by Prange (cf. again Ref. [32^1. 

V. CONCLUSIONS 

We showed by a numerical evaluation of the Lindstedt 
series up to order n — 700 that a previously assumed 
violation of Greene's criterion was ungrounded. The as- 
sumption that allowed the restricted series (|15|) was fal- 
sified for orders n > 200. The resulting critical con- 
stant did not correspond to the SM, but is still true for 
the SSM. From our numerics, we conclude that, for the 
golden mean, the SSM critical constant is strictly higher 
than the SM. This seems to be generally true for all wind- 
ing numbers, but is specifically difficult to proof for the 
golden mean where both constants are very close. Still, 
the numerics till order n — 700 do not give a complete 
convergence. An evaluation that would compare to the 
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FIG. 3: Singularities in K of the function h(x, K) for the SM 
without small divisors for x varying in [0,1]. The radius of 
convergence p corresponds to the value x = 1/4, while the 
critical constant K c = 1 corresponds to x — 1/2. The curve 
is symmetric with respect to both the real and the imaginary 
axes. 



accuracy of Greene's method would rely on an prohibitive 
computational effort. 

In addition to this analysis, we have purposed a new 
model that appears when the small denominators in the 
SM and SSM are suppressed. We show that this model 
maintains many features of the SM and SSM. However, 
it has an analytical solution and corresponds to Kepler's 
equation in case of the SM. Also here, the analogue of 
the SM has a lower value than the one of the SSM. More- 
over, surprisingly the radii of convergence are lower than 
the true models for golden mean winding numbers. This 
proofs that the golden mean winding numbers are re- 
markably resistant to the small denominator effect and 
falsifies a common misconception that the small denom- 
inator problem is the dominant and only mechanism for 
the breaking analyticity transition. The fact, that the 
simplified model still has a transition with a value even 
lower than the true maps for the golden mean, shows that 
this is not the case. 

Finally, we studied the full analyticity domain for the 
two models. Also here, there are striking differences be- 
tween the simplified- and the true maps. Similar to the 
true maps, the set of singularities form a natural bound- 
ary. However, whereas the SM and SSM all the singu- 
larities in K for the function g(K,x) are present for any 
value of x, the situation is different for the simplified 
maps. The simplified SSM has only one singularity in 
the complex iGplane for the function h(K, x) at each 
value of x. The simplified SM has for each value of x 
two singularities symmetric with respect to the real axis, 
except for the singularities on the real axis for x = and 
x = 1/2. The closed natural boundary is retained after 
gathering all singularities for all x. 

This natural boundary is perfect circle in case of the 
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simplified SSM, while it is more stretched curve for the 
simplified SM with a discontinuity on the real axis at 
K = ±1. This shows that the radius of convergence of 
the simplified SSM equals its critical constant, like was 
found for the true maps with golden mean winding num- 
bers. In contrast, the simplified SM has a critical con- 
stant of K c = 1 that is higher than its radius of conver- 
gence. In that respect, the simplified SM resembles more 
the true SM with winding numbers close to rational val- 
ues. Also this is a bit of a surprise, as one would expect 
the contrary, but is consistent with the trend mentioned 
above. Somehow, it is almost as if the model, in which 
all small divisors were eliminated, still suffers more from 
this effect that the SM with the golden mean. 

Therefore, we believe that the study of these kind of 
simplified analytical models are a worthy prerequisite for 
the understanding of the SM, SSM and FK models and, 
in particular, the influence of the small denominators. 
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APPENDIX A: INDUCTION PROOF 



The second relation of Eq. (|19fl is more difficult. One can 
write 



n—m-\-l 



n \ (n — n' — 1 



to * — ' \n' I \ to — 2 

n' = l 

'(ra'-l)/„ _ M(n-n'-m+l) 



n— m+1 



- T -r 

77 Z — ✓ 7 7/ I 



(n-n' - 1)! 



to ^-f n'\(n — n')\ (to — 2)!(n — n' — to + 1)! 
x n'^'- 1 \n~n'f n - n '- m+1 ^ 

n — rn+1 



- Y - 



m ' n'\(m — 2)!(n — n' — m + 1)! 

71—1 

n-m+l . , 

1 ^-^ n\ In — m+l 

to ^— ' (to — 2)\(n — to + 1)! V nf 

n' — l N 

x n / (»'- 1 )(n-n , ) Cn ~ n '~ m) 

to — 1 ^-^ / n \ I n — m + 1 
to ^— ' \ m — 1 / \ n' 

n — 1 v 7 v 

x n y '(n — n p ' 

m — 1 / n \ \ / n — to + 

TO \ TO — 1 / ^— ' \ 7l' + 1 

v 7 n'=0 V 

x (n' + 1)™' (n - n' - i)(»-"'-i-"0 

(to — l)(n — TO + 1) / T7, 



1. Proof of Eq. (J20j) 



TO \m — 1/ * — ' \ n 

X (n' + l)"' _1 (n - ri - l){n-n'-X-m) 



(A3) 



Assuming that relation (121 H is true up to some Taylor 
order n— 1, then the first relation of Eq. (|19fl for n yields: 



Using Abel's identity [34 



i?(n, 1) = nV 



(n-1) 



n — 1 — m 



(n-2)! 



^ (m — 1)! (n — 1 — to)! 



(n-l) m (n-2)! 



' (n — 2 — to)! to! 

n - 2\ 
\ m J 

m=0 v 7 

((n-lj + lf-^n"- 1 , (Al) 



= " £("- !) 



(x + y)(x + y - ahf 1 = ^ ( fc Jxy(a: - a/c)^ 1 



fc=0 



x [y-a(fi-*)] n , (A4) 



with fe = n', n = n — m, a = — 1, x = l,y = m — 1 yields 



where in third equality we have used the binomial theo- 
rem 



3=0 



(A2) 



n'=0 ^ 



(A5) 



/ 1 i \n — m—n — 1 



x (1 + n')" _1 (n-n'- 1) 



11 



Hence 



R(n, m) 



(to — l)(n — to + 1) / n 



7 n'=0 



n — to 
n' 



= (n — m + 1) I | // 

\m — 1 

n — to + 1 



n—m—l 



n (to — l)!(n — to + 1)! 



(n-1)! 



(to — l)!(n — to)! 



n- 1 

TO — 1 



(A6) 



APPENDIX B: TREE FORMALISM 



indeed they are the only labels of the trees, and their val- 
ues uniquely determine the elements of T n ,n- Therefore 
we can rewrite XV: as 



2ni 2™ ^ 

m\-\- . . .-\-m n =n—l 

l (-l)™^"- 1 



l 



TOl! . . . TO„! 



2?ri 2™ n! ' 
where we have used the multinomial theorem 

^ n\ 

m-i ! 



TOi! . . . TO„! 



(B4) 



X^ 1 . . . X n n — [x\ ~\~ • • • ~1~ %n) s 



which extends the binomial theorem to n > 2; see 
§24.1.3. 



1. Proof of Eq. <f20l» 



Proof of Eq. lf25|) 



First of all, by defining a — 2ttx and u(a) — 2ng(x), 
one can write the functional relation that the function 
u(a) has to satisfy as u(a) + {K/2i) exp(ia + iu{a)) = 
for the SSM and u(a) + K sin(a + u(a j) = for the SM. 
Note that in the case of the SSM the function u, which 
in principle depends on two parameters K and a, is in 
fact a function of the only parameter rj = Ke la . 

In terms of the function u(a) the functional equation 
(gj becomes, for the SM, 

2u(a)-u(a+2wl)-u(a-2Trl) = -K sin(a+u(a)), (Bl) 

in which we recognize Eq. (1.4) of Ref. [2(j, with e = K 
and lo = I. For the SSM we have the same equation with 
the sine function replaced with (2i) exp(ia + iu{a)). 
Then we can envisage the same tree expansion as in Ref. 
[20|: see formula (2.2), where, to make a relation with 
the notations we are using now, k and v are what we are 
denoting with n and k, respectively. Moreover 7 (nui v ) = 
—oj vi , hence it is —1 in our case, and one has v v = 1 for 
the SSM and v u e {±1} for the SM. At the end we find 



1 (~l) r 
2m 2 n 



i)er„ 



Val(tf), Val(tf) = TT ——, v 



In the case of the SM, without small divisors, we can 
still use formula (|B2|) . but now one can have v u = ±1. 



1 (—l)n+{n-k)/2 
■y-n _ \ J 

k ~ 2ni 2" 



(n-k)/2 



v 



in 1 



y 

f-^ TOi!...TO„! 
rn 1 + . . . + m n — n — 1 

I (_-\\n-\-(n-k) fc n_1 



(B6) 



2th 



((n-fc)/2)!((n + fc)/2)!' 



m„+l 



for k = n — 2p, with p = 1, . . . , n. 

By using the definition Eq. (|16f) of radius of conver- 
gence, we see that we have to compute the maximum 
over k of \X%\. 

As k = J2uet> v u we see that, first, k can assume only 
the values —n, —n + 2, —n + 4, . . . ,n — 4, n — 2,n (so that, 
in particular, (n ± fc)/2 is even), and, second, in order 
to have a contribution to XV: we have to put (n — k)/2 
mode labels v u equal to —1 and the remaining (n + k)/2 
mode labels equal to 1. Moreover for any tree d £ T n ,k 
we can write 



u6i? 



m„ + l 



(B2) 

where the trees the branching numbers m u and the set 
of trees T n k of order n (that is with n nodes) and with 
momentum k flowing through the root line (that is such 
that X^ugiJ Vu = ^) are defined as in Ref. |20| . 
In the case of the SSM, Eq. (|B2|I reduces to 



ue-9 



uGi9 



n 

(B7) 

which inserted into Eq. (|B2I) gives, by using again the 
multinomial theorem, 



xrn 

-A- — 



i (-i) T 

2ni 2" 



E 



Val(tf), Val(tf) = TT -^-r, 



1 (_l)n+(n-fe)/2 n 

~ ~ ' (n - k) j2 



2iri 



mGiJ 



E 



(B3) 

as v u = 1, and the sum over trees of order n can be writ- 
ten as a sum over all possible configurations of branching 
numbers {m u }ueD with the constraint J2 u ea m " = n ~ 1- 
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